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1. Introduction 



The idea that deterministic theories and quantum theories may be very closely related, 
has been advocated by this author for a number of years[l]. To understand this relation, 
some issues needed to be further clarified. It has always been obvious that determinis- 
tic, time reversible systems can be cast in a quantum mechanical formalism, where the 
evolution is described by a hamiltonian that in no fundamental way differs from other, 
obviously quantum mechanical hamiltonians. However, in many cases, different choices 
of the hamiltonian describe the same physical model, in particular if the time variable is 
restricted to integral time values (in some units). On the other hand, one would like to 
implement locality of a deterministic system, such as a cellular automaton, by postulating 
a hamiltonian that can be regarded as the integral (or sum) of a local hamilton density. 
This latter procedure is not straightforward. We found interesting examples where this 
can be done, but a systematic formalism still seems to be lacking. 

Recently, we elaborated a scheme that enables us to map systems with pairs of integers 
onto other systems where these pairs are replaced by pairs of real but non commuting 
numbers[2]. This mapping requires the machinery of quantum mechanics. States of the 
discrete system may well be postulated to evolve completely deterministically, but this 
is no impediment against treating them as states in a Hilbert space. After our mapping 
is applied, much of their deterministic nature becomes hidden, and the system takes the 
shape of a quantum theory. Indeed, it is a real quantum theory, including Born's inter- 
pretation of the squares of the amplitudes as representing probabilities. Consequently, 
we arrive at a theory of quantum wave functions (which is what quantum states look like 
when defined on the real line), evolving completely in line with an ordinary Schrodinger 
equation, yet it is a deterministic theory in disguise. The importance of this mapping is 
its potential usefulness in quantum field theories. There, we have hamiltonians that are 
integrals of hamilton densities. We set out to follow the following strategy: turn a simple, 
exactly soluble cellular automaton into an exactly soluble quantum field theory, derive its 
hamiltonian, and subsequently add interaction terms. The hamiltonian will receive small, 
local corrections so that now a cellular automaton can be mapped onto an interacting 
quantum field theory. From there, find our way to the Standard Model. 

We claim that, in spite of its quantum mechanical appearance, the Standard Model 
may also be a deterministic system in disguise. The author is aware of the fact that most 
of his readers will not be prepared to jump so such a conclusion; in particular those who 
have Bell's inequaities[3][4] in mind will not be inclined to accept the idea. But then we 
would invite the reader to study the transformation presented in this paper. Here, we 
display how the mapping goes for non-interacting, massless bosons in one space- and one 
time dimension. The resulting quantum field theory is restricted to a discrete space-time 
lattice, but the lattice length, a , may be taken to be as small as we wish, and indeed 
imposing a lattice cut-off is standard practice in quantum field theory. In more realistic 
theories, such a lattice is suspected to replace some much more sophisticated ultraviolet 
cut-off induced by the gravitational force at the Planck scale. 

The quantum hamiltonian of the theory is modified by the lattice artifacts in a way to 
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be specified later; in any case, there will be left-moving particles and right-moving particles 
that show no interactions at all. Being restricted to wave functions on the lattice, we see 
right-away that these particles represent quantised units of information moving to the left 
and to the right, both with the speed of light. 

Our cellular automaton behaves just like this. It is defined by having two sets of 
integers: Q(x,t) are defined on the same lattice sites, while P(x,t) are defined on x 
integer but t is integer plus \ . These Q and P variables will now also be seen to send 
discrete packages of data to the left and to the right, just like our bosonic model, and 
for this reason one could already suspect that a mapping between the automaton and the 
bosonic model should exist. This is the mapping that will be constructed in this paper. 
The mapping is far from trivial, as the reader will surely discover. 

Initially, we intended also to describe the introduction of mass terms and interaction 
terms, along the lines sketched above, but in the 1+1 dimensional system we encountered 
an obstacle against this: both the cellular automaton and the quantum field theory feature 
a Goldstone symmetry[5] forcing the excitations to be massless. Breaking the Goldstone 
symmetry would require non-local interactions of a type we would rather avoid, while it is 
known that interacting quantum field theories in 1+1 dimensions cannot have Goldstone 
symmetries [6]. We conclude that the introduction of interactions and masses will require 
more advanced techniques or generalization to more dimensions, first. 

Perhaps this is not necessary. Our mapping may have more direct implications for 
string and superstring theories. These do contain massless bosons on a 1+1 dimensional 
world sheet, so our observation may imply that determinism and discreteness may be 
considered in a systematic way for string theory. This could lead to new models, or at 
least new interpretations of the existing models. 



2. The continuum quantum field theory 

Consider first the standard, massless, bosonic field theory in 1+1 dimensions, described 
by a lagrangian and a hamiltonian, 

£ = -\{d t q 2 - d x q 2 ) ; H = J &x{\p 2 + \8 x q 2 ) , (2.1) 

where we use the symbol p(x) to denote the canonical momentum field associated to the 
scalar field q(x) , which here obeys p(x) = d t q(x) . The equal-time commutation rules 
are usually: 

[q(x),q(y)\ = \p(x),p(y)\ = ; [q(x),p(y)\ = i5(x - y) . (2.2) 

If we now include the time variable by writing q(x, t) and p(x, t) in Heisenberg notation, 
the solution of the field equations can be written as follows: 

a L (x,t) = p(x,t) + d x q(x,t) = a L (x + t); (2.3) 
a R (x,t) = p(x,t) - d x q(x,t) = a R (x-t). (2.4) 
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In terms of these variables, the hamiltonian is 



H= J dx\{{a L {x)f + {a R {x)f) , (2.5) 
and using k to indicate momentum variables, we have in momentum space 

H = J™ dk\ ((a L \-k)a L (-k) + a R \k)a R (k)^) , (2.6) 

k,k' >0: [a L (-k), a L (-k')\ — , [a L (—k), a L \—k')} = 2k5(k — k!) . (2.7) 
[a R (k), a R (k')} = , [a R (k), a R \k')} = 2k5(k - k!) . (2.8) 

In our notation, a L,R (k) are the annihilation and creation operators apart from a factor 
\/2k , so the hamiltonian (2.5) can be written as 

POO 

H= dk (kN L (—k) + kN R (k)) . (2.9) 
Jo 

where N L ' R (^fk)dk are the occupation numbers counting the left and right moving par- 
ticles (we subtracted the vacuum values for the energy). The energies of these particles 
are equal to the absolute values of their momentum. All of this is completely standard 
and can be found in all the text books. 



3. Notation 

In the previous Section, we still used (more or less) conventional notation, since everything 
was in the continuum. On the lattice, however, we will pick a somewhat less standard 
convention, as it will simplify the calculations there significantly. In the sequel, we will 
use 

capital Latin letters, N, P, Q, X, • • • , to indicate integer valued fields, (3.1) 
lower case Latin letters, a, p, q, ■ ■ ■ , to indicate real numbers, and (3.2) 

lower case Greek letters, a, rj, £, A, • • • , for fractional numbers, (3.3) 

the latter being usually confined to the interval (— |, |] . An exception is the lattice 
coordinates x, t , which are in lower case because they are not used as field variables. 

Frequent use will be made of the number 

e = e 2n w 535.49 • • • , so that e 2nia = e ia , and e a = l if Z G Z . (3.4) 

If we have integers Q, P, ■ ■ ■ , we will often associate a Hilbert space of states to these: 
| Q, P, • • •) . Then, there will be operators i]q, rjp, • • • , defined by 

e iNr lQ\Q,P, ...) = \Q + N,P, ■■■) ; e iNr lr\Q,P,...) = \Q,P + N,-.-) . (3.5) 
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In general, these operators will have eigenvalues restricted to be in the interval (— |, |] . 
The matrix elements of the r/g operators can readily be calculated [2]: 

We will need to know the commutator between the operators r/q and <3 : 

(Qilfog, QMQ2} = (Q 2 - Qi)(Qil^lQ 2 ) = ^ (<W - (-1)*-*) . (3.7) 

Because we wish to keep factors 2n in the exponents (to be absorbed when we use 
e instead of e), our commutation rules for p and q operators (which are real numbers) 
will be 

[q, p\ = i/2n , (3.8) 
which means that it is h rather than h that we normalize to one. 



4. Quantum field theory on the 1+1 dimensional lattice 

Inserting a lattice cut-off for the UV divergences in quantum field theories is also standard 
practice. Restricting ourselves to integral values of the x coordinate, we replace the 
commutation rules (2.2) by 

[q(x),q(y)\ = [p+ (x) , p + (y)] = ; [q(x),p + (y)] = j-J x , y (4.1) 

(the reason for the superscript + will be explained later, Eqs. (4.14) and (A.l)). The 
exact form of the hamiltonian on the lattice depends on how we wish to deal with the 
lattice artifacts. The choices made below might seem somewhat artificial or special, but it 
can be verified that most alternative choices one can think of can be transformed to these 
by simple latticle field transformations, so not much generality is lost. It is important 
however that we wish to keep the expression (2.9) for the hamiltonian: also on the lattice, 
we wish to keep the same dispersion law as in the continuum, so that all excitations must 
move left or right exactly with the same speed of light. 

The lattice expression for the left- and right movers will be 

a L (x + t) = p + (x,t) + q(x,t) - q(x - l,t) ; (4.2) 
a R (x-t) = p + (x,t) + q(x,t)-q(x + l,t) . (4.3) 

They obey the commutation rules 

[a L , a R }=0; [a L (x), a L (y)} = ±± if y = x ± 1 ; else ; (4.4) 
[a R (x), a R (y)} = if y = x ± 1 ; else . (4.5) 
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In momentum space, writing 

1 

a L ' R {x) = j\ &Ka L ' R (K)t lRX , (4.6) 

the commutation rules are 

[a L (/t!), a L (-K 2 )} = - k 2 ) sin(27TKi) ; (4.7) 

[a fl (/ci), a R (-/t 2 )] = J <J(k! - k 2 ) sin(27TKi) . (4.8) 

If we want the hamiltonian to take the form (2.9), then, in terms of the creation and 
annihilation operators (4.7) and (4. 8), the hamiltonian must be 

i 

H= [ 2 dK . ™ J a L (K)a L (-K) + a R (K)a R (-K)) . (4.9) 
J sin(27r/c) V / 

Since, in momentum space, Eqs. (4.2) and (4.3) take the form 

a L {n) = + (1 - e- iK )q{n) , a R (/t) = p+(«) + (1 - e iK )q(K) , (4.10) 

after some shuffling, we find the hamiltonian 
l 

H = tt [ 2 dn( K \p + {k)\ 2 + AKt&mtK \q(n) + ± P + (k)\ 2 ) , (4.11) 
J q Vtan7r/^ " / 

where |p + (/t)| 2 stands for p + (/t) p + (—/t) . Since the field redefinition g(a;) + ip + (a;) — >■ q(x) 
does not affect the commutation rules, and 

Urn — ^ = 1 , 4sin 2 (™)|g( K )| 2 |(^)(«)| 2 , (4.12) 

we see that the continuum limit (2.1), (2.6) is obtained when the lattice length scale is 
sent to zero. 

We now notice that the operators a L (x,t) = a L (x + t) and a R (x,t) = a R (x — t) move 
exactly one position after one unit time step. Therefore, 

a L (x,l) = a L (x + 1,0) , a R (x,l) = a R (x - 1,0) , etc. (4.13) 

and now we can use this to eliminate p + (x,t) and q(x,t) from these equations. Writing 

p + (x,t)=p(x,t + ±) , (4.14) 

one arrives at the equations 

q(x,t + l) = q{x,t) + p(x,t + \) ; (4.15) 

p(x,t + l) = p{x,t-\) + q(x-l,t)-2q(x,t) + q(x + l,t) . (4.16) 
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We now see why we had to shift the field q(x, t) by half the field momentum in Eq. (4.11): 
it puts the field at the same position t + \ as the momentum variable p + (x, t) . 

Thus, we end up with a quantum field theory where not only space but also time is 
on a lattice. The momentum values p(x, t + 1) can be viewed as variables on the timelike 
links of the lattice. 

At small values of k, the hamiltonian (4.9), (4.11) closely approaches that of the 
continuum theory, and so it obeys locality conditions there. For this reason, the model 
would be interesting indeed, if this is what the cellular automaton is going to match. 
However, there is a problem with it. At values of k approaching k — >■ ±| , the kernels 
diverge. Suppose we would like to write the expression (4.9) in position space as 

H = J2 M \s\ (a L (x) a L (x + s) + a R (x) a R (x + s)) , (4.17) 



X, s 



then M s would be obtained by fourier transforming the coefficient ttk/ sin(27r/t) on the 
unit interval for k . One obtains 



M = — e— =±r° gj ~ Ek= ° ^ ^ S = 6Ven (4 18) 

s 7_i +A sin(2^) * \ log(2A) + E J-i)/2 i if s = odd 

where A is a tiny cut-off parameter. The divergent part is 

^( io §x) Y,(- i y x ~ v { aL ^ aL ^+ aR ^ aR ^) = 

^(iog^)((E(- 1 )^ L ^)) 2 + (E(- 1 )^)) 2 ) • ( 4 - 19 ) 



Also the kernel 4Ktan7Tfi; in Eq. (4.11) diverges as k — > ±| . Keeping the divergence 
would make the hamiltonian non-local, as Eq. (4.19) shows. We can't just argue that the 
largest k values require infinite energies to excite them because they do not; according to 
Eq. (2.9), the energies of excitations at momentum k are merely proportional to k itself. 
We therefore propose to make a smooth cut-off, replacing the divergent kernels such as 
4«;tan7r/t by expressions such as 

(4/ttan7r/s:)(l-e A2( |- K)2 ) , (4.20) 
where A can be taken to be arbitrarily large but not infinite. 



5. The cellular automaton 

Our cellular automaton is a model defined on a square lattice with one space dimension 
x and one time coordinate t , where both x and t are restricted to be integers. The 
variables are two sets of integers, one set being integer numbers Q(x,t) defined on the 
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lattice sites, and the other being defined on the links connecting the point (x, t) with 
(x, t + 1) . These will be called P(x, t + \) , but they may sometimes be indicated as 

P + (x, t) = P-(x, t + 1) = P(x, t + \) . (5.1) 

The automaton obeys the following time evolution laws: 

Q(x,t + 1) = Q(x,t) + P(x,t+l) ; (5.2) 
P(x,t+\) = P{x,t-\) + Q(x-l,t)-2Q(x,t) + Q(x + l,t) , (5.3) 



just analogously to Eqs. (4.15) and (4.16). It is also a discrete version of a classical field 

4< 



theory where Q(x, t) are the field variables and P(x, t) = %Q{x, t) are the classical field 



momenta. 

Imagining the x, t lattice to be very large, one might be interested in studying the 
statistical properties of the Q and P fields. To this end, we introduce Hilbert space, just 

as a tool. The basis elements of this Hilbert space are the states {Q(x, 0)}, {P + (x, 0)}^ . 
If we consider a superposition of such states, we will simply define the squares of the am- 
plitudes to represent the probabilities. The total probability is the length- squared of the 
vector, which will usually be taken to be one. At this stage, superpositions mean nothing 
more than this, and it is obvious that any chosen superposition, whose total length is one, 
may represent a reasonable set of probabilities. The basis elements all evolve in terms of a 
permutation operator that permutes the basis elements in accordance with the evolution 
equations (5.2) and (5.3). As a matrix in Hilbert space, this permutation operator only 
contains ones and zeros, and it is trivial to ascertain that statistical distributions, written 
as "quantum" superpositions, evolve with the same evolution matrix. 

Indeed, in what follows later, we will make frequent use of operators in this Hilbert 
space, such as the operators f]Q{x) and i]p(x) , which are defined exactly as in Eq. (3.5) 
in Section 3, but now at each point x at time t = (the operator r]p(x) acts on the 
integer P + (x, 0) ). 

The time variable t is an integer, so what our evolution equations generate is an 
operator U{t) obeying U{t\ +t 2 ) = U{ti)U{t2) , but only for integral time. In principle, 
it would be easy to find an operator H , to be called "hamiltonian" , such that 



U(t) = e~ iHt ; H = J2 ^tr^ U ^) - f/ ("^)) • ( 5 - 4 ) 



71=1 

This equation is obtained by Fourier analysis: when — \ < x < \ , we have 



x 



00 fl 2(— l) n_1 

a n sin 2-kux ; a n = 2 / i xsin27ra; = — . (5.5) 

71=1 2 



The problem with this hamiltonian is that 

1. It is not unique: one may add any integer to any of its eigenvalues; and 
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2. It is not extensive: if two parts of a system are spacelike separated, we would like 
the hamiltonian to be the sum of the two separate hamiltonians, but then it will 
quickly take values more than 1 / 2 , whereas, by construction, the hamiltonian (5.4) 
will obey \H\ < l / 2 ■ 

Thus, by adding appropriate multiples of integers to its eigenvalues, we would like to 
transform our hamiltonian into an extensive one. The question is how to do this. 

Indeed, this is the question that forced us to do the investigations described in this 
paper; the hamiltonian of the quantum field theory considered here is an extensive one, 
and also naturally bounded from below. 

At first sight, however, the similarity between this automaton and the quantum field 
theory of Section (4) may seem to be superficial at best. Quantum physicists will insist 
that the quantum theory is fundamentally different. 

Our procedure will force us first to compare the left-movers and the right-movers in 
both theories. In Appendix A, some of the preperties of systems obeying the evolution 
equations (5.2) and (5.3) are listed. Here also, it is found that the operators a L (x) 
and a R (x) have the desirable properties that they are left movers and right movers, and 
furthermore, they can directly be expressed in terms of the original field variables. In the 
automaton, these are integers, so, the results of Appendix A imply that we must write 
them as: 

A L (x + t) = P + (x,t) +Q(x,t) -Q(x-l,t) ; (5.6) 
A R (x-t) = P + (x,t) + Q(x,t) -Q(x + l,t) . (5.7) 

Note, that these are fields that transport streams of data to the left or to the right, but 
at this point the resemblance to elementary particles may still seem to be remote. 



6. The mapping 

We now claim that all states of the cellular automaton of Section 5 can be mapped onto 
the quantum field theory states of Section 2, with a few exceptions: the mapping is nearly 
though not quite one-to-one. Several ways to describe this mapping were tried, and we 
found the most transparent one to be as follows. We first consider the operators f]p(x) 
and t)q(x) , as described in Sections 3 and 5. They will be identified with operators in 
the quantum field theory, by observing the following. 

If we have operators and Pi obeying the commutation rules 

h, Pj) = £ s a ' ( 6 - x ) 

then we can readily derive the commutation properties of the exponentiated operators: 

£ iapi e ibqj = e iab5ij ^ibqj £ iapi _ /g_ 2 ) 
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This means that, if a and b are both integers, e mPt and e %hqi commute. If we now define 
the operators rjQi and r] P j as follows, 

e iv Qi = e -i Pi ; e iv Pj = e ni ; ( 6 3 ) 

while furthermore restricting r/Qi and r/pj all to lie within the interval ( — \ \ ] , then they 
are uniquely defined, and they all commute. Thus, we can define all operators t]q(x) and 
r)p(x) of our cellular automaton as (minus) the fractional parts of the field momentum 
operators p + (x, 0) and the field operators q(x, 0) of our quantum field theory. Note, that 
e ~t P +(x) ra j ses the field q(x) by one unit and e iq ^ raises the field p + (x) by one unit, so 
indeed the operators that raise or lower fields q(x, 0) by one unit, will now do the same 
for the variable Q(x,0) of the automaton, and the operators the raise or lower p + (x, 0) 
in the field theory, do the same for the variables P + (x, 0) in the cellular automaton. In 
the automaton, these raising and lowering operators all commute, but so they do in the 
quantum field theory. 

Note, that the identification of the raising and lowering operators in the automaton 
with the ones in the quantum field theory, is restricted to be valid for the states at time 
t = . To find how the identification is done at other times, we use the Heisenberg 
notation: operators are time dependent, while states are time independent. 

What has been achieved so-far is that if we have identified one state in the cellular 
automaton with one state in the quantum field theory, then we can identify all other 
states, simply by applying the raising and lowering operators as often as we want 1 . How 
do we identify one state? We can rephrase this question: 

We identified the fractional parts of the fields q(x, t) and p + (x, t) with cellular automaton 
operators. How do we identify the integral parts? 

We can't just identify them with the variables Q(x,t) and P(x,t) because the time 
evolution laws will round off the fractional parts differently. And therefore, this turns out 
to be not so easy. Take the cellular automaton "vacuum" state 

\n) = \{Q(x,0) = 0, P + (x,0) = 0}) . (6.4) 

Could this be mapped onto the ground state of the quantum field theory's hamiltonian? 
This is probably incorrect since the state \Q) is orthogonal to all states obtained by 
raising or lowering one of its Q or P + values, but when applying one of the operators 
e iNp+(xfl) or e iNq(xfl) ^ energy ground state of the quantum field theory, these in 
general do not annihilate this state. 

Could \VL) be mapped onto the state Y\ x ipo{q{x ) Q)) , where ip is the wave function 
derived in ref [2], which is orthogonal to all those states that would be obtained by 
raising or lowering Q and/or P + values? Again, this would probably be wrong, because 
we cannot prove that this state is time-independent (even when we restrict ourselves to 
integer time values), whereas the state \Q) is. 

1 Indeed, this part of the identification is valid in many more cellular automaton models that can be 
identified with quantum field theories, also in higher dimensions; it is the next part of the argument 
however, that we do not know exactly how to generalize. 
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The correct mapping is obtained by viewing Hilbert space as the product of the space 
spanned by the left-movers and that spanned by the right-movers. We do the mapping 
for the left-movers and the mapping for the right-movers separately. 

So let us first concentrate on the left movers. The cellular automaton is then described 
by the integers A L (x + t) , as defined in Eq. (5.6). In that Hilbert space, we also have the 
operators tja(x + t) , defined by 

e iN ^\{A(y)}) = \{A'(y) = A(y) + N5 xy }) , (6.5) 

where the superscript L was temporarily omitted. As in Eq. (3.6), we can write the 
matrix elements of t\a '■ 

{A^x^a^^x)) = —{5 Al { x ),A 2 {x) - l ) a 2 (x) - A^x) ' ^ 6 ' 6 ^ 

while r]A(x) commutes with all operators A(y) and i]A(y) when y^x. 

From Eq. (6.6), we derive the commutator between tja and A . This is as in Eq. (3.7): 

(A^rjA, A]\A 2 ) = JL (6 Ai a 2 - (-1)*"*) . (6.7) 
This may be written as 

[ VA , A] = ±(i- IVaXVuI) ; (A\iPa) = (~1) A ■ (6.8) 

1^) is a very special state that is defined at every spacelike point x . So we have 

[ VA (x), A{y)\ = £5 xy (l - \Mx))(Mx)\) ■ (6-9) 

We call \iPa(x)) an edge state, since it lives on the edge of the interval in r\ space: 
Va(x) = ±| . 

Our next attempt to construct the mapping consisted of omitting these bothersome 
edge states. Then, 

[ VA (x),A(y)]l^5 xy . (6.10) 

This would allow us construct operators a L (x) that obey the commutation rules of the 
quantum field theory, Eq. (4.4). Re- inserting the superscript L : 

a L (x) = A L (x) + i] L A [x - 1) . (6.11) 

and similarly we coud define the right-movers. The good point about this attempt is, that 
the time evolution just shifts these operators one step to the left or to the right, without 
mixing them. Since in the cellular automaton, A L [x) and r\\ (x) shift the same way, this 
mapping commutes with the time evolution. 

Eq. (6.11) is not demonstrably false, since the constraint that we limit ourselves to 
the states that are orthogonal to the edge states \^a(x)) does commute with the time 
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evolution. However, it is an unnatural constraint. The edge states would have finite 
energy, in general, so using these operators a L and a R as creation and annihilation 
operators might lead to wrong results. 

A superior mapping procedure is the one introduced in Ref. [2]. We know how the 
r] operators are transformed, but to express the operators P + and Q themselves, or 
equivalently, the operators A L (x) and a R (x) , we have to establish a phase function in 
i] space, which is now an infinite-dimensional torus, that is, the product of all circles 
defined by the periodic operators r)p(x,0) and t]q(x,0) . The procedure goes exactly as 
in Ref. [2]. 



We replace (6.11) by 



d d 



where <f(ff^) is a phase function that depends on all values of f]\{x) , but it is the same 
function for all points x . One easily derives that, indeed the commutation relation (4.4) 
holds — as long as tp is non-singular, and the term rj^(x — 1) should not jump. 

But, of course, r}\{x — 1) does jump, when we restrict it to stay within the interval 
(— |, \ \. This is what we need the phase (p for. It is postulated to obey the periodicity 
conditions 

<p({r&(x) + S X , X1 }) = ¥>({ftf (*)}) + itf (a* + 1) . (6.13) 

If, in Eq. (6.12), we substitute x — 1 = x\ , we see that this expression is now strictly 
periodic in 77 space, so now we expect no contribution from edge states. 

Yet there is a contribution from an edge state, and this is because the periodicity 
condition (6.13) cannot hold everywhere, since there is a clash. If we go around the torus 
in two different directions, to return to the same point, a full phase rotation results. Thus 
9 7 ({ ? ?i( ;r )}) nas vortices trapped in it, and these cause new singularities. 

The function is easy to construct [2] . Write 



00 



^}) = E^a(^ + i),^)); rfr, J2 c -2<*-0-**», (6 . 14) 

X K = — OD 

where r and •p are both real functions of rj and £ . This is a special case of the elliptic 
theta function -$ 3 , and it can also be written as a product [7] instead of a sum 2 . Due to 
this property, one can also write the function (p(r),£) as follows: 

pfaS) = ( ar g(! + €^-^-2) + arg(l + e-^-^-2)) . (6.15) 

K=0 

There is a duality property: 

<p(v,Z) + <p(Z,v) = vt- (6-16) 



2 The function r(j7,£) here differs by a factor e 2^ from the function used in Ref [1]. 
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Eq. (6.12) simplifies into 



« L (*) = -£^y + g^(v(v L A^ + l),V L A^))-^V^-l),V L A ^))) • (6-17) 

We can now compute [2] the matrix elements of a L (x) when sandwitched between two 
cellular automaton states |{A L (:r)}) : 



{A[\a L (x)\A^) = A L (x)5 + 5 

- 5 



(_ 1 jAAL( x )+AAL(x+l) +HAA L( x+1 j 

x, x +i A#(x) 2 + A#(i+1) 2 

(_ 1 )AA^(x)+AA^(x-l)+l iAA L( a . _ X ) 

x,x-i A#(x) 2 + AA L (x-l) 2 



(6.18) 



where 5 is the identity matrix in A L space (so that Af(x) = A%(x) = A L (x) ), while 
S\y jZ stands for the matrix that keeps all values of A\ (x) equal to Af (x) , except at the 
points x = y and x — z , where we define 

AA L (y) = A L 2 (y) - Af(y) , AA L (z) == Afc) - Af(z) ; (6.19) 

(the diagonal element of S\ VjZ , where the terms in (6.18) would be ill-defined, are taken 
to be zero). 

It is now important to check the commutation rules for the field a L (x) . Since the 
operators A L (x) refer to the states of the automaton, they all commute with one another, 
and since the additional terms in (6.18) are diagonal in rj space, they also commute with 
one another (this is trivially seen to be true in Eq. (6.18) because the matrix elements 
only depend on AA L , not on the A\ 2 (x) themselves). The non-commuting parts come 
from the cross terms. Since a matrix element of [X, A L (a;)] in A L space is the matrix 
element of X multiplied by AA L (x) , one finds 



[a L (x), a\x+l)} = £(l - \^(x,x + l))(^(x,x + 1)|) , (6.20) 

where the new edge state is defined by 

(A\x), A L (y)\^ L A (x,y)) = (-l^M+^M . (6.21) 

This edge state originates from the flux singularity in the phase function tp(r),£) located 
at the corners r] = ±|, £ = ±| , and it can be replaced but not removed. We claim that 
this state is far less harmful than the edge states encountered previously. The point is 
that any state that is not orthogonal to any of these edge states, carries infinite energy. 
This we see by calculating the expectation value of the hamiltonian (4.17). Plugging the 
matrix elements (6.18) of a L (x) in here, we see a logarithmic divergence of positive terms 
(when s = 0), and negative terms (when s = 2), but the positive ones dominate since 
M > M 2 . The odd s terms do not diverge. 
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We conclude that the edge state contribution in Eq. (6.20) can be neglected whenever 
states with finite energy are considered. In our construction, all eigen states of the op- 
erators a L (x) are unique. Repeating the same procedure for the a R operators gives us 
a complete mapping of all finite energy states of the quantum field theory onto states of 
the cellular automaton. The mapping commutes with the time evolution, since both the 
a L (x) operators and the a R (x) operators, just like the ones of the cellular automaton, 
merely shift to the left or to the right as time evolves. An important restriction is that 
we should only look at those elements of the automaton states that are orthogonal to the 
edge states: 



where (P,Q) = ( A L (x), A L (x + 1) ) or (P,Q) = ( A R (x), A R (x + 1) ) for any point x 



on the lattice. 

Curiously, this condition excludes the pure basis elements of the automaton such as 
the state \Vt) , which is why in a more sophisticated theory such a constraint might be 
undesirable. We observe that the exceptional states, the edge states (6.21), consist of a 
single state at every site. As such, it may describe fermions on the lattice. In a next 
publication, it will be explained that the bothersome constraint is lifted if we map the 
cellular automaton of Section 5 onto a super symmetric 1+1 dimensional quantum model, 
with one bosonic field together with one fermion. 

In the quantum field theory, we found another mild restriction: all states with particles 
in the very highest energy mode, k — ±| , contribute non-locally to the hamiltonian, since 
the kernel (4.9) diverges there. This is why we propose the cut-off (4.20). 

We emphasize once more: the important reason why this mapping is of interest is, that 
now we know the hamiltonian, Eq. (4.17), where we can plug in the operators a L ' R (x) 
since we know their matrix elements (6.18). 

This result is quite non-trivial. Since this hamiltonian correctly transports the left 
movers to the left and the right movers to the right, by one step in every unit of time, 
also the integral parts A L (x) and A R (x) move correctly, as well as their fractional parts, 
T)\{x — 1) and rj R (x + 1) . So this is the correct hamiltonian that generates the evolution 
law of the automaton. Being identical to a quantum field theory hamiltonian, it no longer 
suffers from any positivity problem, which was an important difficulty in previous work[l]. 

7. Interactions 

Up till this point, our construction might be looked upon as a mere mathematical curiosity, 
but physically rather boring, since both the cellular automaton and the quantum field 
theory decribe nothing but non-interacting particles. The left-movers do not react at all 
on the presence of right-movers and vice versa. 

In principle, introducing interactions at both sides of the mapping seems to be straight- 
forward. Suppose that we make a minor modification in the cellular automaton's evolution 




(6.22) 
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law. After every move as the time clock counts, it is checked whether some rather im- 
probable combination of cellular field variables occurs, for example: at any value of x , 

Q(x) = 0, P + (x) = 5, Q(x + 1) = 17, Q(x - 1) = 3 . (7.1) 

Only if this situation occurs, P + (x) is augmented by one unit. Another condition is 
considered where P + (x) decreases by one unit. We write this action as a permutation 
matrix: 

V = e~ iHint , H int = \J2 W n \x) . (7.2) 

X 

We already had the hamiltonian H = XL^o(^) , the hamiltonian (4.17). To write the 
effect of both interactions in a hamiltonian form, we use Baker-Campbell-Hausdorff[8]: 

e A e B = e A+B + llA,B] + - ) (73) 

where the dots represent an infinite series of commutators. The fact that these are com- 
mutators is very important, because it means that the net result can again be written as 
the integral over a hamiltonian density, and this means that locality of the theory is not 
dramatically affected. Now, we took the condition (7.1) to be sufficiently uncommon so 
that its effect may be expected to be small. Therefore, perturbation expansion in powers 
of the parameter A in Eq. (7.2) makes sense; the higher terms only contribute of the con- 
dition is met in several ways, which will reduce the numerical values of those contributions 
even more. We do see now that the quantity B + \[A, B] + ■ • • in Eq. (7.3) generates 
an 'interaction hamiltonian' in the quantum field theory. It is important to observe that, 
with this modification, the hamiltonian not only represents a theory with interaction, but 
it still can be written as the sum over hamiltonian densities. Since the hamiltonian den- 
sity is corrected by finite correction terms, it continues to have well-defined lower bounds: 
the vacuum state. 

However, in this particular case, there is a serious complication that practically inval- 
idates this procedure: to express the original fields q(x,t) of the cellular automaton in 
terms of the a L and a R operators that we constructed, one has to invert Eq. (4.10) , and 
for small k values this blows up. Clearly, only the derivatives of the q fields occur in the 
original model, not the q fields themselves. So, in our special case, it seems that we can 
only handle interactions where the interaction hamiltonian is a function of dq/dx rather 
than q itself. However, if we try to get interactions this way, we cannot lift the Goldstone 
symmetry[5] that made the particles massless. Not only might our original intention to 
create interaction and mass terms be difficult to realize in this model, so that any at- 
tempts to create mass terms might give rise to unacceptable non-localities; unfortunately, 
this also jeopardises our attempts to get any interactions at all, because interacting field 
theories in 1+1 dimensions are known not to allow for a Goldstone realization of a contin- 
uous symmetry [6]. We believe that this difficulty is an artifact of the model considered; 
perhaps more sophisticated models can be found, in particular in higher space dimensions, 
where such difficulties do not arise. As for this paper, we do not expand any further on 
the topic of interactions. 

Our main intention with this paper was to demonstrate in what way cellular automata 
can actually behave as quantum field theories and vice versa. 
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8. Discussion 



It is difficult to explain more clearly the theory that there may be determinism behind 
what is called quantum mechanics today, than by pointing to the mapping displayed in 
this paper. Every operator in the quantised field theory has its counterpart in the variables 
describing a cellular automaton. What exactly this means concerning the interpretation of 
the EPR paradox[4] and Bell's inequalities [3], is difficult to phrase precisely. An attempt 
was made by the author in Ref. [9]. For one thing, our theory is definitely not a hidden 
variable theory of the type that Bell considered in his paper [3]; however these arguments 
clearly failed to convince many people. Perhaps the best argument consists of stating that 
our theory is nothing but an extrapolation of the mathematics phrased here, claiming that 
there is no logical obstacle against continuation of these models towards systems that more 
and more resemble the Standard Model of the Elementary Particles. 

Yet there are a few subtleties. A subset of the cellular automaton states, being all 
states that are not orthogonal to the edge states, have to be left out of the mapping, 
as they would have infinite energy in the quantised field theory. This constraint is not 
insignificant, because it means that for all left moving modes and for all right moving 
modes, at every pair of points, Eq. (6.22) implies that we must postulate that 



and the cellular ground state, \Q) = |0, 0, • • •) , does not obey this constraint. We must 
conclude that either the quantum field theoretical description of cellular automata only 
applies to automata that are in rather highly excited modes, far away from the trivial 
state , or we have to add fermions to the bosonic theory. 

Secondly, the resulting quantum field theory has a hamiltonian that approaches the 
continuum theory smoothly as the lattice length a tends to zero, but the lattice edge 
states, the particles at maximal momentum, k — ±| , contribute to the hamiltonian non- 
locally, so that the kernel M s of the hamiltonian is non-local after all. We would insist 
that this feature would be unnoticeable for observers in a universe who have no access to 
particles with such a high momentum. 

Note however that we do have all states of the quantised field theory in our cellular 
automaton, as we reproduce Eq. (2.9), describing all possible occupation numbers of the 
particles. Note also that the cellular automaton model itself is entirely local, so if we use 
that as our "ontological" underlying theory, the apparent signals of some non-locality in 
the effective bosonic theory cannot be used to explain what some observers regard as a 
necessary non-locality in hidden variables theories for quantum mechanics. 

In principle, interactions can be introduced, by adding small modifications to the 
cellular evolution rule, as explained in Section 7. However, it is also explained there that 
there are pathologies in 1+1 dimensions. These are caused by the Goldstone symmetry 
that we built in. As is well-known [6], interacting quantum field theories cannot realise 
Goldstone symmetries in 1+1 dimensions, while we also found it difficult to break the 
symmetry. This also obstructs the generation of mass terms. 






A L (x),A L (x+l) 
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Our mapping may have interesting implications for string theory. Bosonic strings are 
harmonic theories on a one plus one dimensional string world sheet. If we introduce a 
lattice in the world sheet, strings become small solid 'rods' stitched together, but in our 
model also target space, the 10 or 26 dimensional space-time, was made discrete. This 
means that our cellular automaton is relevant for string theory, describing discrete strings 
on a space-time lattice. [10] 

A. Left and right movers 

The cellular automaton as well as the quantised field theory will be postulated to obey 
the same equations on the lattice. In both systems, we locate the field variables q(x, t) 
on the regulare lattice sites (x, t) where both x and t are restricted to be integers. The 
momentum variables are located conveniently on half-odd time coordinates, but sometimes 
it will be more practical to locate them on integer spots. Therefore, we write 

p ± (x,t) =p(x,t±\) . (A.l) 

While in the quantized field theory these variables are operators with a continuous spec- 
trum of eigenvalues, obeying the usual commutation rules (4.1), the cellular automaton 
is described by replacing q(x,t) — > Q(x,t) and p(x,t) — > P(x,t) , where Q and P only 
take integer values and all commute. In this Appendix we keep the notation q and p . 

The evolution equations are Eqs. (4.15) and (4.16). It is sometimes useful to have the 
complete solution of these equations forward in time: 

t 2t-l 

9(0, f) = ^p + (2n-l-t, 0) + ^(-ir^g(m-t, 0) ; (A.2) 

n=l m=l 
2t 

p + (0,t) = J2(-l) m (p + (m-t, 0) + 2q(m-t, 0)) - q(t, 0) - q(-t, 0) . (A.3) 

m=0 

The q and p fields can be split up in left-movers and right-movers: 

q(x,t) = q L {x + t) + q R {x-t) ; p(x, t + \) = p L (x + t + \) + p R (x - t - \) , (A.4) 

which obey the equations 

p L (x + \) = q L (x + 1) -q L (x) ; p R (x - \) = q R (x - 1) - q R (x) . (A.5) 

However, we will often consider the more fundamental left and right movers a L , a R , 
defined by inverting the equations (A.4) : 

a L (x) = p L (x + \)+p L (x-\) = p + (x,0) +q(x,0) - q(x - 1,0) ; (A.6) 
a R (x) = p R {x + \)+p R {x-\) = p + (x,0) + q(x,0)-q(x + l,0) . (A.7) 

It is easy to verify from the basic equations (4.15) and (4.16) that the rhs of equations (A.6) 
and (A.7) will roll to the left and to the right, respectively. In momentum space, a L and 
a R will be associated to the particle creation and annihilation operators of the field theory. 
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Assuming the commutation relations (4.1), the operators a L and a R will obey easy 
commulation rules themselves, as given in Eqs. (4.4) and (4.5), while the commutation 
rules for q L , q R , p L and p R will be non-local. This is why it would be difficult to intro- 
duce interactions that are expressed directly in terms of the q(x) and p + (x) operators. 
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